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Abstract 

, We show that the Maxwell equations for arbitrary inhomogeneous me- 

' dia are equivalent to a single quaternionic equation which can be con- 

| sidered as a generalization of the Vekua equation for generalized analytic 

""-■■■« . functions. 

■ 

A| 1 Introduction 

-t— > ■ 

S Quaternionic reformulation of the Maxwell equations in a vacuum is quite well 
known (see, for example, Q, jfj], 0, Q). The system of Maxwell equations 

> '. 

.!_, , div D = div B = 

a) 

rot E = -d t B, rot H = 9 t D, 

where D = EqE and B = ^qH (eq and /zo are the permittivity and permeability 
of free space), is equivalent to the quaternionic equation 

(-d t + iD)7=0, (2) 
Co 

where Co is the speed of light in a vacuum, D is the Moisil-Theodoresco operator 
(see the definition below) and f is a purely vectorial complex quaternion. The 
equivalence between (Q) and (||) can be established by putting 

Re / = y/e~oE and Im / = ^JIqH. 

The quaternionic approach to Maxwell's equations for homogeneous media was 
intensively used in a number of works (e.g., H, H, pOj, p3[), but even the 
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question as to how to write the Maxwell equations for arbitrary inhomogeneous 
media in a compact quaternionic form remained open until recently (an attempt 
in this direction can be found in |5|, Section 4.5]). In such a reformulation 
was proposed in the case of a time-harmonic electromagnetic field and in [ p"l"[ 
for the time-dependent case. Here we make one additional step which leads us 
to Maxwell's system in the form of a single quaternionic equation, which can 
be considered as a generalization of the well known in complex analysis Vekua 
equation describing generalized analytic functions |^5| . 



2 Preliminaries 

We will consider continuously differentiable functions of four variables (t, xi, X2, x$) 
with values in the algebra of complex quaternions H(C). By D we denote the 
operator D = J>2k=i ^kdk- Here dk — g§- and are the imaginary quaternionic 
units. Let us notice the following property of the operator D. Let tp be a scalar 
complex function and g be an H(C)-valued function. Then 

D[ V >-g] = D[y]-g + ip-D[g\. (3) 

Taking into account that D[<p] = grad</? = iid±ip + iidi^p + i^d^ip and assuming 
that ip is different from zero we can rewrite ([|) in the form 

iD + E**)g=*D[ V .g]. (4) 

We will use the following notations for the operators of multiplication from the 
left-hand side and from the right-hand side 

a Mg := a • g and M a g := g ■ a, 

where a € H(C). The usual complex conjugation we denote by Vectors 
from C 3 are identified with purely vectorial complex quaternions. Note that 
for an H(C)-valued function g — go + if the action of the operator D can be 
represented as follows 

Dg = — div if + grad go + r °t if ■ 



3 Maxwell equations 

We assume that the relative permittivity e r and the relative permeability fi r of 
the material arc differentiable functions of coordinates e r = e r (x\, x%, X3) and 
fi r = fi r (xi, X2, X3). The permittivity and the permeability of the medium are 
introduced as follows 

£ = So£ r and ^ — /ioMr- 
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Then Maxwell's equations for an inhomogeneous medium have the form 

rotH = £0 t E+j, (5) 



rotE = -/z<%H, (6) 
div(eE) = p, (7) 

div (jiH) = 0, (8) 

where all the magnitudes are real. Equations (0) and (||) can be written as 
follows 

„ grade „ p 
divE+ < ,E >= - 

e e 



and 



divH+<i^,H> =0 , 



where < •,• > denotes the usual scalar product. Combining these equations 
with (|^) and (0) we obtain the Maxwell system in the form 

^ E=< grad £E> _p 

e e 

and 

DH=< gradM H >+ea t E+j. (10) 
A* 

Let us make a simple observation: the scalar product of two vectors ~j? and if 
can be represented as follows 

< y, Y >= M + m*)-?. 

Using this fact, from (^) and ([!(]) we obtain the pair of equations 

(D + i^£) B = _ 1 M ^E - p dt H - I (11) 



and 



Note that 



(D+ lg^a> )H= _l M ^ H + £atE+ . (12) 



1 grad e grad ^/i 

2 e = Jl ' 
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Then using (||), equation (|Tl]) can be rewritten in the following form 



-^£(0 • E) + E • ^ = -/i<9 f H--, (13) 

J £ £ 



where 

gradv^e 



Analogously, (|l^) takes the form 

• H) + H • 7? = ea t E+j, (14) 

v M 

where 

Introducing the notations 

t := y/iE, Jt := vT^H 

and multiplying ( jl3| ) by y^i and (O) by y//I we arrive at the equations 

( J D + Af 7 ')£ > = -i9^--^, (15) 
c V £ 

and 

(£> + M-*)Tt = -d t ~t + Vw, (16) 

where c = 1/^/eJl is the speed of propagation of electromagnetic waves in the 
medium. 

Equations @ and @ can be rewritten even in a more elegant form. 
Consider the function 

~? + 0t 

Let us apply to it the quaternionic Maxwell operator 

1 



d t + iD. 

c 

We obtain 

(-dt + iD)~T = ~d t ~? - DTt + i(-d t Jt + D~t). 
c c c 

For the real part of this expression we use equation ( |l6| ) and for the imaginary 
part equation (p~5]). Then we have 



Note that 



Hence 



{-8t + %Dj~? = -i(M^t + iM^lt) -JJii-lL, (17) 

c y £ 



~? = \{T + 7*) and 7? = I(7-7«). 
2 2i 



Let us notice that 

. , grade _. , gradW 

e + u = and e — u — — - — , 

c W 

where W — \J fi/e is the intrinsic wave impedance of the medium. Denote 

_> := grad^c ^ w ;= gradVg 

vc Vw 



Then 

M^~t + iM^rt = ~(M^y+M w y*). 

From ((T^) we obtain the Maxwell equations for an inhomogeneous medium in 
the following form 

(-d t + iD)J-M^7-M iW 7* = -(Vmj + ±L) (18) 
c V e 

(compare with (||)). This equation is completely equivalent to the Maxwell 
system (^)-(^) and represents Maxwell's equation for inhomogeneous media in 
a quaternionic form. 



Remark 1 Equation \ldj can be considered as a generalization of the well 
known in complex analysis Vekua equation describing generalized analytic func- 
tions Jlq/. Recently in JZ^ / using the L. Bers approach J^/, J^/ another quater- 
nionic generalization of the Vekua equation was considered. Probably some of 
the interesting results discussed in JT^ / can be obtained for ( fJ^ j also. Then their 
physical meaning would be of a great interest. 
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